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ABSTRACT 

We have compiled a large sample of 151 high redshift (z = 0.5 — 4) galaxies selected at 24^m 
{S 2 4 > 100 fjjy) in the GOODS-N and ECDFS fields for which we have deep Spitzer IRS spec- 
troscopy, allowing us to decompose the mid-infrared spectrum into contributions from star formation 
and activity in the galactic nuclei. In addition, we have a wealth of photometric data from Spitzer 
IRAC/MIPS and Herschel PACS/SPIRE. We explore how effective different infrared color combina- 
tions are at separating our mid-IR spectroscopically determined active galactic nuclei from our star 
forming galaxies. We look in depth at existing IRAC color diagnostics, and we explore new color- 
color diagnostics combining mid-IR, far-IR, and near-IR photometry, since these combinations provide 
the most detail about the shape of a source's IR spectrum. An added benefit of using a color that 
combines far-IR and mid-IR photometry is that it is indicative of the power source driving the IR 
luminosity. For our data set, the optimal color selections are S^so/S^ vs - S8/S3.6 and Sioo/S^ vs. 
•58/53.6; both diagnostics have ~10% contamination rate in the regions occupied primarily by star 
forming galaxies and active galactic nuclei, respectively. Based on the low contamination rate, these 
two new IR color-color diagnostics are ideal for estimating both the mid-IR power source of a galaxy 
when spectroscopy is unavailable and the dominant power source contributing to the IR luminosity. 
In the absence of far-IR data, we present color diagnostics using the WISE mid-IR bands which can 
efficiently select out high z (z ~ 2) star forming galaxies. 



1. INTRODUCTION 

In the current narrative of galaxy evolution, star for- 
mation and the growth of supermassive black holes are 
intertwined. The star formation ra te density of the Uni- 
verse peaks from z ~ 1 — 3 (e.g.. iBouwens et al.l l2009t 
iMagnelli etahl 120111: iMurphv et all 1201 1[) . an epoch in 
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which the black holes within the center of ma ssive galax- 
ies are simultaneously b uilding up their mass (IWall et al.l 
I2005t iKellv et al.l l2010K There is compelling evidence 
that the growth of black holes and the buildup of stellar 
mass is linked, though the processes that regulate, and 
ultimately quench, the simultaneous growth of the stellar 
mass and black hole are not yet fully disentagled (e.g., 
IMullanev etld1[20Dl . 

To study the properties of active galactic nuclei (AGN) 
and star- forming galaxies at z ~ 1 — 3, when the star 
formation rate in the most massive galaxies begins to 
decline, it is necessary to first identify systems likely 
harboring an AGN. Because AGN are much more lu- 
minous in the X-ray than star-forming galaxies, X-ray 
detection pro vides one of the best m eans of identifying 
an AGN (e.g., Alexander et al. 2003). Deep field surveys 
with the Chandra X-ray Observato ry have exposed an 
AGN population out to z ~ 5 (e.g.. IBrandt et al.|[200ll: 
iGiacconi et al.ll2002T ). However, detailed X-ray spectral 
analysis shows that t he majority of sources are obscured 
by gas and dust (see IBrandt fc Hasingerl [20051 for a re- 
view), and there is likely a sizeable fraction of AGN that 
remain undetected by Chandra surveys, as evidenced by 
the fact that nearly half of the X-ray bac kground is un- 
resolved at > 6 keV (jWorslev et alj [20051 ) . and the ratio 
of observed obscured to unobscured AGN at high red- 
shifts is lower than what is found for c omparably lumi- 
nous AGN in the local Universe (e.g., iTreister fc Urrvl 
I2005T ). 

Due to the incompleteness of surveys conducted at 
X-ray wavelengths, we must employ alternate methods 
insensitive to dust obscuration to identify the presence 
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of an AGN. The infrared portion of the spectral energy 
distribution (SED) shows clear signs of both AGN and 
star formation (SF) activity. Based on the presence of 
polycyclic aromatic hydrocarbons and continuum ther- 
mal dust emission, mid-infrared spectra can be decom- 
posed into the relative contributions of SF and AGN 
activity (e.g . . lLaurent et all l2000t [Armus et all 120071; 
Sauna et all 120071; iPope et al.l T2008; Ki rkpatrick et al.l 



2012aD . At 1.6 fim, SF dominated galaxies will ex- 



hibit a stellar bump due to emission from older stellar 
populations, whereas as the dust surrounding luminous 
AGN heats up, it will radiate into the near-IR and mid- 
IR, creating a pure power-law spectrum at these wave- 
lengths. Finally, the peak and shape of the far-IR emis- 
sion depends on the temperature of the dust, which be- 
comes warmer as the AGN grows m ore luminous (e.g., 
lHaas et al.ll200l iSanders et al.l[l98l . Due to AGN sig- 
natures from the near-IR to the far-IR, infrared color- 
selection techniques are a promising way to select out 
AGN missed by X-ray surveys, or when X-ray data are 
unavailable. 

In this work, we explore the IR color space for a 
sample of 151 high redshift luminous infrared galaxies 
(LIRGs, Lm = 10 11 — 1O 12 L ) and ultra luminous in- 
frared galaxies (ULIRGs, Ljr > 10 12 L Q ) at high redshift 
(z ~ 0.5 — 3.5) with deep Spitzer mid-IR spectroscopy 
and a suite of multiwavelength photometry spanning 
3.6 — 500 /im. Several studies have explored u sing Spitzer 
IRAC colors as a means to separate AGN (|Lacv et al.l 
l200i IStern et al.ll20"ol IDonlev et al1l20T2l) . For the first 
time, we explore these various diagnostics using a sam- 
ple of mid-IR spectroscopically determined AGN and SF 
galaxies. In light of the emerging Wide- held Infrared Sur- 
vey Explorer (WISE) photometry, we also discuss how 
well mid-IR diagnostics separate our spectroscopic AGN 
and SF sources. The focus of the paper, however, lies 
in combining near-, mid-, and far-IR photometry to dis- 
tinguish our AGN from the SF galaxies. The far-IR is 
tracing the bulk of the IR luminosity, and if a galaxy's 
mid-IR power source is also affecting the far-IR, we ex- 
pect that combining both portions of the spectrum will 
be a useful diagnostic. 

We build on a previous paper (jKirkpatrick et al1l2012aL 
hereafter Paper I) in which we separate our sample into 
AGN and SF dominated, based on mid-IR spectral de- 
composition, and analyze the average IR SEDs of each, 
including mid-IR spectral features, IR luminosity, and 
dust temperatures. We have spectroscopically deter- 
mined the nature of the mid-IR power source for our 
individual sources, and we have created SEDs that rep- 
resent the average features of our SF and AGN sources. 
Now we look for advantageous color-color cuts that can 
be used to select out AGN candidates when spectroscopy 
is unavailable. We explore different combinations of col- 
ors using both the photometry of our individual sources, 
as well as redshift tracks from the composite SEDs, to 
determine color combinations that separate AGN from 
SF galaxies. 

The paper is laid out as follows: in Section 2, we de- 
scribe details of our sample, the mid-IR spectral decom- 
position of individual sources, and the composite spectral 
energy distributions (SEDs) we create from this sample. 
In Section 3, we explore the efhcacy of existing IRAC 



color selection techniques for high redshift sources. In 
Section 4, we present new color diagnostics based on com- 
bining photometry from the Spitzer Space Telescope and 
the Herschel Space Observatory. We apply our new di- 
agnostics to the broader GOODS-N and ECDFS fields in 
Section 5, and we end with our conclusions in Section 6. 
Throughout this paper, we assume a standard cosmology 
with H = TOkms^Mpc -1 , fl M = 0.3 and fl A = 0.7. 

2. DATA 

A full description of our sample and the composite 
SEDs we create from it is given in Paper I. Here, we 
summarize the main details and results. 

2.1. Multiwavelength data 

Our sample consists of 151 high redshift galaxies from 
the Great Observatories Origins Deep Survey North 
(GOODS-N) and Extended Chandra Deep Field Survey 
(ECDFS) holds. We include all sources in these fields 
that were observed with the Spitzer IRS. While this sam- 
ple contains a diverse range of sources depending on the 
goals of each individual observing program, the overlying 
selection criteria is that each source must be detected at 
24 /im with a flux of 5*24 > 100 /xJy, since anything fainter 
will not be observable with the IRS in less than 10 hours. 
More details on this database of IRS sources in GOODS- 
N and ECDFS can be found in our data paper (Pope et 
al. in preparation). 

The GOODS fields have been extensively surveyed 
and are rich in deep multiwavelength data includ- 
ing: ground-based imaging in the near-IR (J and 
K bands) from VL T /ISAAC (iRetzlaff et al.l l2010h an d 
CFHT/WIRCAM (IWang et all 120101: ILin et al.l 12012ft: 
Chandra 2 Ms X -ray observations (|Alexander et aLi r2003: 
ILuo et al.ll2008h : 3.6, 4.5, 5.8, 8.0 ^m from the Infrared 
Array Camera (IRAQ on Sytizen IRS peak-up observa- 
tions at 16 yitm (iTeplitz et al.ll2011l) and MIPS imaging at 
24 and 70 fim (Ma gnelli et al.ll2011ft . Recently, GOODS- 
N and GOODS-S have been surveyed with the GOODS- 
Herschel Open T ime Key Program (P.I. David Elbaz, 
lElbaz et al.1 [2011 using both the PACS and SPIRE 
instruments providing deep photometry at five far-IR 
wavelengths: 100, 160, 250, 350, and 500 /mi. For the 
present study, we combine space-based imaging from 
Spitzer and Herschel to obtain 12 photometric band- 
widths spanning the near-IR to the far-IR. 

For sources lacking a detection at the Herschel wave- 
lengths, we extract a measurement of the flux density 
and associated uncertainty for each source directly from 
our images. The images are in units of mJy/beam, so 
we use the 24 /Ltm prior positions to find the appropriate 
pixels for each galaxy. We do not take a measurement 
when a source looks too blended on the image itself. In 
Paper I, we present complete Herschel photometry for 
our IRS sources, as well as indicating which sources have 
measurements and which have detections at the Herschel 
wavelengths. 

We performed spectral decomposition of the Spitzer 
IRS mid-IR spectrum for each source in order to dis- 
entangle the AGN and SF compo nents. We follow 
the technique outlined in detail in IPope et all (2008) 
which we summarize here. We fit the individual spec- 
tra with a model comprised of three components: (1) 
The SF component is represented by either the local 
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Rest Wavelength (\im) 

Fig. 1. — Composite SEDs (and associated uncertainties) for each 
of our subsamples (see Table[lJ created by stacking photometry and 
spectroscopy in the near- and mid-IR and fitting a two- temperature 
model to the far-IR. The composites have been offset to allow for 
easy comparison. The SF templates lack uncertainties in the range 
~ 18 — 30 (im as we lacked data in this wavelength range, and so we 
interpolate between the mid-IR and far-IR portions of the spectra. 
Composite SEDs are publicly available. 

starburst composite of IBrandl et al.l ((2006) or simply 
the mid-IR spectrum of the p rototypical starburst M82 
(|Forster Schreiber et all 120031) - with the SNR, wave- 
length coverage and spectral resolution of our high red- 
shift spectra both give equally good fits to the SF compo- 
nent of our galaxies. (2) The AGN component is deter- 
mined by fitting a pure power-law with the slope and nor- 
malizatio n as free parameter s. (3) An extinction curve 
from the iDraine et al.l (|2003t ) dust models is applied to 
the AGN component. The extinction curve is not mono- 
tonic in wavelength and contains silicate absorption fea- 
tures, the most notable for our wavelength range being 
at 9.7 ^m. We fit all three components simultaneously 
and integrate under the PAH and continuum compo- 
nents to determine the fraction of the mid-IR luminosity 
5 — 12 /im) from SF and AGN activity, respectively. 
For each source, we quantify the strength of the AGN 
in terms of the percentage of the total mid-IR luminos- 
ity coming from the AGN continuum component. Based 
on this mid-IR spectral decomposition, we find that 38 
(25%) out of our sample of 151 galaxies are dominated 
(> 50% of luminosity) in the mid-IR by an AGN. 

2.2. Galaxy Classifications 

To more thoroughly compare the mid-IR spectral prop- 
erties and full IR SEDs within our sample, we divide our 
sample into four primary subsamples based on the results 
of the mid-IR spectral decomposition. First, each galaxy 
is classified as either SF- or AGN- dominated based on 
having < 50% and > 50% AGN contribution in the mid- 
IR, respectively. We further divide the SF galaxies into 
two bins: z ~ 1 (z < 1.5) and z ~ 2 (z > 1.5). The AGN 
sources are likewise separated into two bins: those with 
measurable 9.7 /im silicate absorption (hereafter referred 
to as silicate AGN) and those without (hereafter referred 
to as featureless AGN). We are unable to further classify 
four AGN sources as they lack spectral coverage in the 
relevant range (9 — 10 /im) to determine whether they 
exhibit silicate absorption. We refer to these as unclas- 
sifiable AGN in the relevant figures. Our four primary 
subsamples are listed in Table [1] along with their median 
rcdshifts and observed frame median 24, 100, and 8 /urn 



flux densities. 

While the majority of our sources that are classified as 
SF-dominated based on the mid-IR spectra have a neg- 
ligible (< 20%) contribution from an AGN, the AGN 
dominated sources exhibit varying degrees of concurrent 
SF activity. The featureless AGN (lacking silicate ab- 
sorption) primarily have a very strong AGN continuum 
accounting for 80 — 100% of the mid-IR emission, whereas 
the silicate AGN have a more uniform distribution of 
AGN fraction (50 - 100%) with some silicate AGN also 
having weak PAH features. 

2.3. Composite SEDs 

To assess the average properties of our primary four 
subsamples of galaxies (due to the small number, we do 
not examine the unclassifiable AGN in detail) , we created 
composite SEDs from 0.3 — 600 /im rest-frame by combin- 
ing data from ground-based near-IR; Spitzer IRAC, IRS, 
and MIPS (24 /mi, 70 /mi); Herschel PACS and SPIRE. 
We reject any source with less than 3 measurements or 
detections in the far-IR bandwidths, or any source with- 
out mid-IR data in the range 6.4 — 7.5 /im, as this is 
the range used to normalize the individual SEDs. Af- 
ter normalization, we stacked flux densities in bin sizes 
of ~ 0.1 /tm below 20 /im. Above 20 /im, we fit all data 
points with a two-temperature modified blackbody. Full 
details of these composites are found in Paper I. The 
composites are shown in Figure [1] with each composite 
offset on the y-axis to allow for easier comparisons. Com- 
posite SEDs are publicly available^- 

The two SF composites are remarkably similar in shape 
and have most emission from cold dust. The feature- 
less AGN composite is nearly a pure power-law until 
~ 20 /tm, and then is relatively flat from ^ 20— 100 /im. 
The silicate AGN SED is a power-law in the near-IR, has 
weak PAH features and silicate absorption at 9.7 /im, has 
warm dust emission around 20 /tm, and has a cold dust 
component peaking at the same wavelengths as the SF 
SEDs. The difference in shapes between the two AGN 
SEDs and the SF SEDs suggests that IR color diagnostics 
could be useful for separating AGN from SF galaxies. 

3. IRAC COLOR-COLOR DIAGNOSTICS 

Some of the most well used IR color diagnos tics for sep- 
aratin g AG N and SF ga l axies a re presented in lLacv et al.l 
p00l and iStern et all pOOll) and utilize IRAC colors. 
The motivation behind an IRAC selection technique is 
that, at these wavelengths, luminous AGN should have a 
monotonically increasing SED, and these power-law col- 
ors will separate AGN from SF galaxies in colorspace. 
With our high redshift photometry and spectroscopy, 
we are able to apply these diagnostics to a large sam- 
ple of mi d-IR spectro s copica lly d etermined AGN an d SF 
galaxies. lLacv et al.l (|2004f ) and IStern et al.l (|2005t) de- 
fine IRAC color-color regions to separate AGN based on 
large surveys of low redshift (z < 0.7) galaxies. As we 
move to higher redshift (z ~ 2), the IRAC bandwidths 
begin to probe the stellar bump; our composite AGN and 
SF SEDs exhibit different shapes in the A < 4 /im region 
of the spectrum, indicating that IRAC color diagnostics 
might also be useful at hi gher redsh ift. We apply the 
diagnostics of IStern et all (j2005[ ) and lLacv etaLl (|2004l ) 

15 http: / / www.astro.umass.edu / ~pope/Kirkpatrick2012 / 
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TABLE 1 

Basic properties of our four sub-samples. 



Subsample 


# of a 
Sources 


Median (Mean) 
Redshift 


(10 12 L ) 


Median S24 c 

(MJy) 


Median Sioo c 
(mjy) 


Median S 8 c 


z ~ 1 SF galaxies 
z ~ 2 SF galaxies 
Silicate AGN 
Featureless AGN 


69 (39) 
44 (30) 
22 (17) 
12 (9) 


1.0 (0.9) [0.8. 1.0] 
1.9 (1.9) [1.8, 2.1] 
1.9 (1.7) [1.6, 2.0] 
1.2 (1.1) [0.6. 1.6] 


0.42 ± 0.17 
2.00 ± 0.71 
1.65 ± 0.54 
0.76 ± 0.07 


370 [260, 570] 
270 [220, 370] 
470 [250, 860] 
1520 [1240, 2300] 


7.9 [4.8, 14.7] 

3.2 [1.6, 5.0] 

5.3 3.0, 10.2] 
9.5 [4.0, 11.7] 


42 [23, 57] 
17 [13, 38] 
64 [21, 140] 
288 [240, 311] 



We list the upper and lower quartile values in brackets next to each calculated median. 

a We list the number of sources in each sub-sample that are used to create the composite SEDs in parent heses . We do not 
include sources with an incomplete IR SED when creating the composites to avoid biasing our results (See M2.3p . 
b Calculate by integrating under the composite SEDs from 8 — 1000 /im. See Paper I for details. 
c Observed frame fluxes. 



z - 1 SF galaxies 




[5.8] - [8.0] (AB) log(S 5B /S 3 . e ) 



Fig. 2. — Color-color diagnostics from Stern ct al. (2005, , left panel) and Lacv ct al. (2004, , middle panel) for selecting AGN dominated 
sources (shaded regions). Our SF galaxies separate nicely with redshift (z < 1.5 are the filled green circles while z > 1.5 are plotted as the 
open green circles; see the online version for color figures). The AGN sources do not separate in these color-color spaces with redshift and 
thus are only plotted according to whether they possess a 9.7 fim absorption feature (open blue squares) or not (filled blue triangles). The 
four AGN we were unable to classify as they lacked spectral coverage at 9.7/im are plotted as the blue crosses. Overplotted are the redshift 
tracks calculated from our composite SEDs for z ~ 2 SF galaxies (purple) and z ~ 1 SF galaxies (red). In the left panel, the z ~ 1 SF SED 
and z ~ 2 SF SED diverge around z ~ 1 due to the difference in the strength of portion of the stellar bump traced by the 3.6 ^m filter. 
As the two AGN composites occupy the same region of colorspace for both the left and middle panels, we only plot the redshift track of 
the featureless AGN composite SED (in pink) on the left and the track for the silicate absorption AGN (orange) in the middle. The right 
panel shows the effective wavelengths of the four IRAC bands at different redshifts on each of our composite SEDs, demonstrating that at 
high redshift these colors sample the stellar bump region. 



to our sample (Figure [5]). Sources in our sample that 
we determined through mid-IR spectral decomposition 
to be dominated by an AGN (> 50%) in the mid-IR are 
plotted in blue (see the online version for color figures), 
and sources dominated by star formation are plotted in 
green. 

In the left panel of Fig. [21 we p lot our sou r ces in 
the IRAC color-color space defined in lStern et al.l (2005) 
with the gray shaded region being the area defined as 
AGN-dominated by the authors. Our SF sources sepa- 
rate cleanly according to redshift and largely avoid the 
gray shaded region. We have also overplotted the redshift 
tracks of our z ~ 1 SF composite SED and our z ~ 2 SF 
composite SED. We calculate the redshift tracks of each 
of our composite SEDs by convolving the SEDs with the 
IRAC bandpass filters (and MIPS, PACS, SPIRE, and 
WISE filters in Q at the appropriate wavelengths for a 
given redshift. The convolution acts to smooth out much 
of the noise in our spectra as the individual filters span 



a fairly large wavelength range. The z ~ 2 SF track con- 
taminates the AGN region around z ~ 1, and in fact, all 
of the z ~ 1 SF galaxies (green filled circles) occupying 
the shaded region have a redshift between 1 and 1.5. The 
divergence of the two SF tracks around z ~ 1 is due to 
the fact that the 3.6 ptm filter is tracing the bluest por- 
tion of the stellar bump, which has differing strengths 
relative to the other IRAC filters for the two SF com- 
posite spectra (see Fig. Q] and the right panel of Fig. [2] 
). Our individual z ~ 1 sources have colors that are lo- 
cated around both the z ~ 1 and z ~ 2 tracks indicating 
a spread in the obscuration of the IRAC colors for these 
sources. The spread of the individual photometry points 
in relation to the composite SED is illustrated in Paper 
I. 

The AGN sources (blue symbols) are not completely 
constrained to the gray shaded region, though our red- 
shift track from the featureless AGN composite SED in- 
dicates they should be on average, illustrating a lack of 
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Fig. 3. — Our sou rces in IRAC color space with the criteria of 
I Donley ct al. (2012) shown as the shaded region. The modified 
criteria do an effective job of selecting power-law sources while 
avoiding contamination from SF galaxies but still miss many of 
our silicate AGN. 

homogeneity in the near- and mid-IR photometry among 
mid-IR classified AGN. Most AGN in our sample do not 
have a clear stellar bump in the near-IR which restricts 
them to the shaded region (only 8% (2) of the AGN 
sources in the grey region have a visible stellar bump) . Of 
the outlying high redshift AGN in this plot, the majority 
(70%) possess a weak visible stellar bump which pushes 
them into the regions of the color space occupied by SF 
galaxies. Eleven (50%) of the silicate AGN, 2 (17%) of 
the featureless AGN, and 1 (25%) of the unclassifiable 
AGN lie outside the gray region. One of the advantages 
of this study is that we are able to investigate how well 
the mid-IR power source determines the near- and far-IR 
power source. Clearly, our mid-IR dominated AGN do 
not uniformly have IRAC colors indicative of an AGN. 
Furthermore, 11 of the AGN contaminating the SF re- 
gion do not have far-IR colors expressive of an AGN (see 
©• 

The middle panel of Fig. [5] is the color space defined 
bv lLacv et al.l (|2004f ) where the shaded region is used to 
identify AGN. The vast majority of our sources occupy 
this region, regardless of their pow er source diagnosed 
from the mid-IR spectrum (see also IDonlev et al.l [20081 
2012). The SF sources again show a clean redshift sepa- 
ration as the IRAC channels sample the stellar bump. As 
neither subsample of AGN similarly exhibits a redshift 
separation, we do not plot them with different symbols 
according to redshift. The redshift tracks of the feature- 
less AGN and silicate AGN composite SEDs occupy the 
same portion of the graph, so we only plot the silicate 
AGN SED track in orange; the track lies in the upper 
portion of the graph. 

We overplot both SF SED redshift tracks (z ~ 1 in 
red and z ~ 2 in purple) as they show an interesting 
separation. At low redshifts (z < 2), the z ~ 1 SF SED 
track lies just outside the shaded region, while the z ~ 2 
SF SED lies inside it, but after z = 2, the tracks lie on top 
of each other. This difference in the color tracks between 
our z ~ 1 and z ~ 2 composite SEDs is likely due to 



the different intrinsic Lir of each subsample. The z ~ 1 
sources are on average LIRGs (£ir ~4x 10 n L Q ) and 
therefore have less dust to obscure the IRAC colors than 
our z ~ 2 ULIRG SED composite; adding more dust to a 
SF galaxy will cause it to shift towards the top- left of this 
plot, which is exactly the shift we see between our z ~ 1 
and z ~ 2 SED tracks. At z > 3, the tracks approach the 
area occupied by the AGN, but as none of our SF sources 
possess a redshift this high, we are unable to determine 
if our sources follow our tracks into the upper portion 
of the graph. Both of the SF SEDs accurately trace the 
redshift separation exhibited by our sources. 

In the right panel of Fig. [5J we show on our composite 
SEDs the effective wavelengths of each of the IRAC band- 
widths at redshifts 1, 2, and 3 (blue, green, and red, re- 
spectively). These bands straddle the stellar bump over 
this redshift range. Both AGN SEDs have a power-law 
shape over these bands causing little evolution in IRAC 
colorspace. 

The high degree of contamination in these IRAC color- 
color diagnostics by hig h redshift SF galaxies moti- 
vated iDonky^eFal] (|2012 ) to cre ate more res t rictiv e cri 
teria than o r iginal ly presented in 
IStern et al.l (|2005D . IDonlev et al 



Lacv et al.l ([2004D and 
()2012[ ) uses a sample 
of AGN identified at optical and X-ray wavelengths to 
construct the new criteria, and now we are able to test 
these criteria using mid-IR spectro scopic AGN. We ap - 
ply the new color cut presented in IDonlev et al.l ()2012D 
to our sample in Figure [3] The authors determined the 
gray shaded region was the most effective at selecting 
out AGN with power-law spectra in the IRAC band- 
widths, and indeed, all of our sources that meet the 
more restrictive criteria have flux densities such that 
S3. 6 < S4.5 < 5*5.8 < Ss.q- Based on our high redshift 
sources, we propose the simple IRAC color cuts of 



log (S5.8/S3.6) > 0.08 
and log (WS4.5) > 0.15 



(1) 
(2) 



(solid lines in Fig. [3]) for selec ting potential power -law 
AGN candidates, similar to the IDonlev et al.l (|2012f ) cri- 
teria. 

The IRAC color selection techn iques, even us i ng th e 
more restrictive cuts presented in IDonlev et al.l (2012), 
still miss a large fraction (39% of the present sam- 
ple) of mid-IR spectroscopically confirmed AGN, specif- 
ically most of the more obscured silicate AGN. Further- 
more, such diagnostics cannot conclusively determine if 
an AGN is significantly contributing to the bolometric 
luminosity of a galaxy. Since an IRAC color diagnostic 
applied at high redshift (z > 1.3) is necessarily based on 
separating sources into AGN- or SF-dominated based on 
the shape of the spectrum in the near-IR, such diagnos- 
tics might not be the most desirable for determining the 
dominate power source of dust obscured galaxies in the 
mid-IR and far-IR regime. 

4. NEW COLOR-COLOR DIAGNOSTICS 

Our large sample of high redshift SF and AGN sources, 
identified with deep mid-IR spectroscopy, and wealth 
of multiwavelength photometry allows us to define new 
color-color diagnostics that are well suited to uncovering 
galaxies harboring an AGN as revealed in the mid-IR 
spectrum, and galaxies with a bolometrically important 
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Fig. 4. — Left - New color-color diagnostic combining submm wavelength Herschel/SPIRE photometry with mid-IR Spitzer data. The 
dark line divides the SF galaxies (green) from the AGN galaxies (blue). We plot all sources according to redshift (filled symbol: z < 1.5; 
open symbol: z > 1.5). We ovcrplot the redshift tracks of our composite SEDs from z ~ 0.5 — 4.0 as this is the redshift range of our sources. 
We show with the dashed grey lines what color thresholds can be used to separate AGN from SF sources if only one color is available. 
Middle - Illustrates where each of the rcdshiftcd photometry points lie on the composite SEDs. Right - Same as the left panel except we 
overplot the SED templates of local galaxies and high redshift SMGs for comparison. 



AGN. We seek to combine multiple portions of the IR 
spectrum that are probing the physical nature of each 
galaxy with different pieces of information. At 3.6 /im, 
SF galaxies at z ~ 1 — 2 will have a stellar bump due to 
the underlying stellar population, while this effect might 
be washed out in luminous AGN, producing a power-law 
shape. Similarly, at these redshifts, the 24 /im filter will 
straddle the PAH complexes at 7.7 and 12.7 /im, which 
will be weakened or absent in a bright AGN. However, 
there is a caveat that at z ~ 1.5, the silicate absorp- 
tion feature present in some AGN spectra, falls into the 
24 /im bandwidth which can produce colors that mimic 
SF galaxies. On the other hand, the bandwidth is large 
enough that we do not expect this to be a significant 
source of contamination. The far-IR should have a dif- 
ferent shape based on the relative amounts of cold and 
warm dust emission present, and SF galaxies will have 
relatively more cold dust than AGN, while AGN have an 
increased amount of warm dust (Paper I). Finally, the 
8/xm filter, at z ~ 1 — 2, covers a relatively featureless 
portion of the spectrum, so when combined with filters 
that are tracing features, should act as a base to distin- 
guish between AGN and SF systems. 

We used available photometry from Herschel and 
Spitzer ranging from 3.6 — 350 /im (observed) and ex- 
plored every possible color combination. In addition, we 
also look at where the composite SEDs lie to get a sense 
of where we should expect to find AGN and SF galax- 
ies, on average. For color-color plots in which the AGN 
were well separated from the SF sources, we calculated 
the contamination rates of each region. We find that a 
color diagnostic spanning the full range of the IR spec- 
trum does the best job of separating both AGN with 
pure power-law spectrum from 1 — 10 /im and AGN with 
silicate absorption from SF galaxies at high redshift. We 
define two new color-color diagnostics which, based on 
having a low contamination and clarity of separation, 
are the optimal diagnostics to employ when a full suite 
of IR photometry exists. 



4.1. S250/S2A vs. S 8 /S 3 . 6 

We find that combining longer wavelength photome- 
try from Herschel with mid-IR photometry from Spitzer 
MIPS/IRAC provides the most reliable separation be- 
tween our mid-IR classified AGN and SF galaxies since 
they probe the widest range of dust properties affected 
by AGN and SF activity. 

Figured]shows S'250/5'24 vs. Ss/S^s- Both colors cause 
a separation of sources with the result that the SF dom- 
inated sources lie in the lower right portion of the graph 
(separated by the diagonal solid line). The SF sources 
separate weakly according to redshift (indicated by open 
and filled symbols) with the lower redshift sources ly- 
ing slightly below the high redshift sources. The redshift 
tracks calculated from the z ~ 1 SF SED and z ~ 2 SF 
SED follow similar evolutionary paths, although again 
they exhibit a separation in the IRAC color, Sg/Sa.e, 
with the z ~ 2 SF track redder than the z ~ 1 SF track 
between z ~ 1 — 2. The SED redshift tracks trace out 
the area occupied by the SF galaxies according to the 
redshift separation exhibited by the sources. 

The AGN are loosely separated by the presence or ab- 
sence of silicate absorption at 9.7 /im in this color-color 
space. We plot both subsamples according to redshift 
(filled symbols are z ~ 1 and open symbols are z ~ 2). 
The featureless AGN do not contaminate the SF region 
and lie further to the left (lower S250/S24) than the sil- 
icate AGN, due to an excess of warm dust emission in 
these sources. We overplot the redshift track computed 
from the featureless AGN SED. The track has some slight 
evolution along the x-axis with redshift, consistent with 
the individual data points. The silicate AGN (filled and 
open squares, according to redshift), on the other hand, 
show no evolution along the x-axis with redshift but show 
some spread along the y-axis. The fact that some sil- 
icate AGN sources contaminate the SF region can be 
attributed to two effects, namely that some individual 
sources possess a weak stellar bump and some sources 
have more relative cold dust emission in the far-IR. The 
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silicate AGN SED redshift track exhibits no evolution in 
either direction and so is plotted as a star. 

Both color axes produce a separation of AGN and SF 
sources. Based on the location of our sources, we calcu- 
late the separation line to be 

log(S 8 /S 3 . 6 ) = 0.74 x log (S 25 o/S 2 4) - 0.78 (3) 

drawn as the bold solid line in Fig. |H In the absence 
of all four wavelengths, we find that the majority of our 
AGN sources satisfy either of the following color criteria 
shown as the dashed lines: 

log(WS 24 ) < 1.30 (4) 
log(S s /S 3 . 6 ) > °- 32 (5) 

The separation along individual axes is particularly use- 
ful if only an upper limit for S250 is available when search- 
ing for AGN. 

We quantify the contamination of the SF and AGN re- 
gions (separated by the diagonal solid line): 10% of the 
sources (11 of 111 total) in the SF region are AGN and 
2 of the 22 sources (9%) in the AGN region are SF dom- 
inated. None of the contaminating AGN are power-law 
dominated. Some of the individual silicate AGN have 
a stellar bump (41%), causing contamination along the 
Ss/Ss.e axis. Futhermore, several of our silicate AGN 
(8, or 42%) have significant amount of cold dust emis- 
sion in the far-IR, producing 1S250 / "5*24 flux density ratios 
greater than 1.3. The high S250/S24 flux density ratio 
correlates with the presence of a stellar bump in the 
silicate AGN, with 32% possessing both. Finally, it is 
worthwhile to note that the contaminating AGN are sig- 
nificantly fainter at 24 /im (median 5*24 = 220 /iJy) than 
the AGN lying in the upper region of the graph (median 
S*24 = 1240 /iJy). Indeed, all but one of the properly 
identified AGN have 6*24 > 300 /iJy. Based on our sam- 
ple of 24 /im bright sources, our diagnostic is optimized 
to select AGN with 5 24 > 300 /iJy, though if the physics 
is similar in lower luminosity AGN, they will have similar 
SEDs and colors, in which case our diagnostics can also 
be used to select out AGN dominated sources. 

S250, in comparison with S 1 ^, is an indicator of the 
relative amount of cold dust in a galaxy. As an AGN be- 
comes more powerful, the relative amount of warm dust 
increases, so that the ratio of decreases. For 

galaxies lacking a significant amount of warm dust, the 
cold dust will be the dominant contributor to the Lir. 
Low S250/S24 ratios in our mid-IR AGN indicate that 
the warm dust, indicative of an AGN, is significantly 
contributing to the far-IR emission, and accordingly, the 
AGN is an important contributor to the bolometric lu- 
minosity. Our diagnostic is therefore more powerful than 
the IRAC diagnostics presented in $3]since it selects AGN 
significantly contributing to the total IR emission. 

Our new diagnostic is a definite im provement o v er the 
IRA C selection c riteria presented in lLacv et all (2004) 
and iStern et ail ((2005), though photometry is scarcer. 
In ou r sample of 151 g alaxie s, 99% have IRAC data so 
both lLacv et all ((2004) and IStern et al.l (|2005f ) can be 
applied, while only 88% have the relevant wavelengths 
to satisfy our new diagn ostic. We verify th e more re- 
strictive IRAC criteria of iDonlev et al.l ((2012:), shown in 
Fig. |3l with our IRS sample, whose mid-IR power source 
has been determined via spectral decomposition. Our 



new diagnostic is a slight i mprovement over th e revised 
IRAC criteria presented in IDonlev et al.l ([2012D . In the 
present case, 35% of our AGN sources are misclassified 
as SF galaxies, while according to the restricted IRAC 
criteria (equations (1) and (2)), 39% are misclassified. In 
addition, with the new criteria, we correctly identify two 
sources as AGN that were misclassified using only the 
IRAC colors. Though the improvements gained in re- 
covering mid-IR AGN by combining colors covering the 
entire IR spectrum is only slight over using IRAC colors 
alone, an added strength of our new criteria lies in the 
fact that mid-IR data from the Wide-field Infrared Sur- 
vey Explorer (WISE) can be easily substituted for Spitzer 
data. The WISE 22 /im and 3.4 /im channels correspond 
to Spitzer 3.6 and 24 /im bandwidths. We have used 
the redshift tracks of our templates to determine that the 
WISE 12 /im channel is the optimal substitute for Spitzer 
8/im (see g2i- 

The middle panel of Fig. @] illustrates the different SED 
tracks shown on the left panel. 250 /im, 24 /im, 8 /im, 
and 3.6 /im at z = 1, 2, 3 are indicated on the composite 
SEDs. S250 traces the far-IR peak of each template while 
S24 traces the mid-IR emission - the different ratios of 
mid-IR to far-IR emission (or warm to cold dust) in AGN 
and SF sources produces the observed color separation 
in S l 25o/'S' 2 4. At the lower wavelengths, Ss/S^.e remains 
fairly constant for both high redshift AGN SED tem- 
plates. As some of our individual AGN sources do show 
signs of weak SF activity as well as AGN activity, sev- 
eral of our AGN possess a noticeable stellar bump which 
causes a spread in Ss/S^.q. Of the 11 AGN sources con- 
taminating the SF region, 60% possess a visible stellar 
bump compared with only 14% of AGN sources in the 
AGN region. The z ~ 1 SF composite SED illustrates 
how a stellar bump would cause a change in Ss/S^.q with 
redshift. 

In the right panel of Fig. 21 we overplot other SED 
templates in our new color-color space, namely the SEDs 
of local galaxies and hi gh re dshift submillimeter galax- 
ies (SMGs. IPope et al.ll2008h . These local SEDs come 
from combining all known data on these sources includ- 
ing available IRAS (12, 25, 60, 100 /im), Spitzer (24, 70, 
160 /im), and SCUBA (85 /im) photometry as well as 
mid-IR spectrosco py (e.g.. iForster Schreiber et al"1l2003l ; 
lArmus et aLll2007l ). Mrk 231 (red) lies in the AGN region 
we defined, although contrary to our high redshift AGN 
SED templates, it does exhibit some evolution along the 
y-axis. M 82 (blue) does not enter the SF region until 
z = 1 which is most likely due to the fact that it's SED 
peaks at a lower IR wavelength than the majority of our 
high redshift sources (see Paper I). Both NGC 6240 (or- 
ange) and the high redshift SMG composite (purple) lie 
in the SF region except at very low redshifts (z < 0.5) 
and follow the same general redshift evolution as our 
sources, that is, increasing S250IS2A color with increas- 
ing redshift. The SMG composite reaches even higher 
5*250 / 5*24 colors than our sample which is consistent with 
their selection at submm wavelengths. The consistency 
of these other local and high redshift templates with our 
new color diagnostics reinforces our confidence in apply- 
ing this color selection to a wider range of IR luminous, 
24 /im bright galaxies at high redshift. 




Fig. 5. — Left - New color-color diagnostic combining far-IR Herschel/PACS photometry with mid-IR Spitzer data. The dark line divides 
the SF galaxies (green) from the AGN galaxies (blue). We plot all sources according to redshift (filled symbol: z < 1.5; open symbol: 
z > 1.5). We also plot the redshifts tracks of our high redshift composite SEDs from z ~ 0.5 — 4.0 which is the redshift range of our sources. 
The silicate AGN composite SED has negligible redshift evolution in this color space and is plotted as the orange star. Middle - Illustrates 
where each of the redshifted photometry points lie on the composite SEDs. Right - Our new color-color plot with templates from local 
galaxies and high redshift SMGs overplotted. 



4.2. Sioo/SW vs. S s /S3.e 

In the absence of longer wavelength SPIRE data, we 
find that we can substitute Sioo for 5*250 &ud we still see 
a nice separation between the SF and AGN galaxies. In 
the left panel of Figure El we use the colors Sioo/SW 
and Ss/Ss.e to define a region that separates our high 
redshift SF and AGN dominated galaxies (diagonal solid 
line). The line of separation is 

log(S 8 /S 3 . 6 ) = 0.208 x log(5 10 o/5 24 ) + 0.105 (6) 

where mid-IR AGN lie above the line. 

We plot the uncertainties on our redshift tracks as the 
hashed lined regions. We opt to not to plot these uncer- 
tainties in the other plots presented in this work (Figs. [2] 
and HJ for the sake of clarity, and the ranges covered by 
the uncertainties in this plot are indicative of the spread 
in the previous figures. As discussed in detail in Paper I, 
the uncertainties on the composites were calculated by a 
bootstrapping technique, which indicates how the scatter 
in the data points affects the calculated median luminos- 
ity by resampling with replacement. The uncertainties 
on our composites are not calculated directly from the 
intrinsic scatter in the data, but are the standard devi- 
ation of the calculated luminosity after resampling the 
data 10,000 times. Therefore, it is not surprising that 
uncertainties on the template tracks do not encompass 
the full spread of all data points, particularly the silicate 
AGN. 

The SF region has only a 12% contamination (11 of 92 
total galaxies) by AGN sources. The SF dominated sys- 
tems show a clear separation with redshift. We overplot 
the redshift track of the z ~ 1 SF composite SED in red, 
and it traces out the evolution exhibited by our individ- 
ual sources. There is only one SF source that lies inside 
the AGN region (causing a contamination rate of 7%), 
and it not only has a high redshift [z — 2.57) but also 
has a 47% AGN contribution to the mid-IR. The silicate 
AGN exhibit no redshift evolution in this color space, 
but the featureless AGN display a weak separation. The 



AGN composite SED tracks (orange and pink) do not 
move much vertically with redshift since these sources 
have a simple power-law shape and lack a stellar bump 
in the rest-frame near-IR. 

The middle panel of Fig. [5] shows three of our com- 
posite SEDs with lines illustrating where the relevant 
photometry bandwidths lie at a given redshift. The 
featureless AGN have a relatively flat spectrum from 
20 — 100 jtzm, so the 100 /im flux does not change with 
redshift whereas the 24 /jm flux decreases, causing the 
evolution along S100/S24. This is not the case for the 
silicate AGN. The slope from the mid-IR to the far-IR is 
relatively constant, producing little change in S100/S24 
at increasing redshifts. 

We have presented a simple IR color-color plot that 
separates our high redshift AGN and SF sources. In the 
right panel of Fig. [5j we test our diagnostic by over- 
plotting the redshift tracks of local templates and a high 
redshift SMG SED. The local AGN Mrk 231 (red dashed 
line) lies well outside the SF region while the starburst 
M 82 (blue) and local ULIRG NGC 6240 (orange) lie in- 
side it for the most part. The SMG redshift track (pur- 
ple) also lies in the SF region as expected since most 
SMG s are star formation dominated (e.g., iPope et "all 
l2008h . 

Substituting Sieo instead of Sioo also works well for 
separating out the AGN (it has the same contamination 
rates mentioned above). However, the SF sources do not 
have a strong redshift separation, although, since S160 
is still probing the cold dust at lower redshifts, there is 
a stronger separation along the x-axis. The S100/S24 
vs. Ss/Ss.q colorspace can be used to select SF galax- 
ies based on redshift as well as looking for AGN. The 
z < 1.5 SF galaxies have higher Sioo/SW ratio since the 
100 /im filter is tracing the cold dust. At redshift of 2, 
the 100 /im filter is now tracing the warm dust, and the 
24 /im flux density is boosted by the 8 /im PAH complex, 
producing a lower Sioo/SW ratio, similar to what is seen 
for the AGN. The similar S W0 /S 2 i colors for both z 2 SF 
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TABLE 2 

Reliability of our new color-color diagnositics. 



Color 
Diagnostic 


# of galaxies 


% (#) Contamination 


SF AGN 


SF region AGN region 


S250/S24, v S%/Sz.q 
Sioo/S'24 V 58/53.6 


102 31 
83 24 


10 (11) 9 (2) 
12 (11) 7 (1) 
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Fig. 6. — We used the WISE transmission filters to create syn- 
thetic 12 and 24 (im photometry for our sources which we com- 
bine with Spitzer photometry of comparable wavelengths to create 
WISE colors. We overplot the redshift tracks of z ~ 1 SF SED 
and silicate AGN SED. The z ~ 2 SF SED and featureless AGN 
SED occupy the same regions as the plotted redshift tracks, so for 
cl arity we do not sho w them. We plot the AGN selection criterion 
of IStern ct al. (2012) as the grey dashed line. Based on the loca- 
tion of our tracks and synthetic photometry, we conclude that the 
primary strength of this mid-IR diagnostic is in selecting z ~ 2 SF 
galaxies. 

galaxies and AGN keeps this color alone from being an 
accurate indicator of the bolometrically important power 
source, unlike <S , 25o/<5 , 24- 

4.3. Far-IR Color Selection 

As an AGN becomes luminous enough to dominate 
the mid-IR spectrum, it can heat the dust in the host 
galaxy causing a shift in the SED to warmer average dust 
temperatures and an increased importance of warm dust 
to the bolometric luminosity. Based on this, we might 
expect that just the Herschel PACS and SPIRE colors 
can can be used to preferen tially select AGN sources. 
lHatziminaoglou et al.l (|2010D searched for a separation 
using S'350/5'250 and £500/5350 and found that in the 
SPIRE bands, their sample of AGN were indistinguish- 
able from the non- active star forming galaxies. As many 
of our sources are not detected or blended at 500 /mi, 
we instead combine S350/S250 with Sieo/S±oo and also 
find that it does not separate the mid-IR classified AGN 
and SF sources. Fig. Q] illustrates that at rest frame 
wavelengths greater than 40 /mi, the far-IR portion of 
the silicate AGN SED and both SF SEDs are all broadly 
consistent in shape explaining the lack of spread in col- 
ors. 

The less pronounced differences in far-IR portions of 
the composite SEDs is also reflected by the failure of 



S250 thresholds alone to preferentially select AGN from 
the larger GOODS survey of galaxies (discussed below; 
see Tableland . Our AGN sources, particularly the 
featureless AGN, are significantly brighter than our SF 
sources at 24 /tm and 8 /mi. Therefore, it is not surprising 
that we found the largest separation between the AGN 
and SF sources when combining mid-IR with far-IR pho- 
tometry. We caution that with only far-IR information 
it is difficult to determine the impact of the AGN on the 
full IR SED. The most reliable selection of AGN candi- 
dates come from combining data from 3.6 /mi to 250 /tm 
( §4.11 §4.2[) . which provides the most detail about the 
shape of IR spectrum. 

4.4. Mid-IR Color Selection 

There is now an abundance of mid-IR photometry 
as a result of two prominent mid-IR space telescopes: 
Spitzer and the Wide-field Infrared Survey Explorer 
(WISE). Past studies have explored using Spitzer color 
combinations to s eparate out mid-IR selected AGN an d 
SF systems (e.g., Ilvison et al.l 120041 : iPope et all 120081 ). 
and emerging studies show that WI SE colors can ef- 
fectively separate IR luminous AGN flE jscnhard t et al.l 
l20ll [Sternet^[2rJll lYan et all 1201 2ft . We are moti- 
vated by these studies to investigate how well the WISE 
photometry can separate our high redshift sources, par- 
ticularly the SF galaxies. WISE has four transmission 
filters centered at 3.4, 4.6, 12, and 22 /tm, and though 
WISE photometry is less sensitive than Spitzer data, it 
has the advantage that it is an all-sky survey and can be 
used to search for high-redshift objects in regions of the 
sky not previously well-studied. 

We do not have WISE photometry for our sources, but 
we can use the WISE transmission filters to create syn- 
thetic photometry using the IRS spectra and the appro- 
priate transmission filters at 12 and 22 /tm. For the 3.4 
and 4.6 /tm filters, we substitute the appropriate IRAC 
photometry. We have applied a small correction to the 
IRAC photometry (0.94 for the 3.6 /im filter and 1.04 for 
the 4.5 /tm filter), which we have calculated by comparing 
the responses of our composite SEDs to each of the WISE 
and IRAC transmission filters. In addition, we calculate 
the redshift tracks of our composite SEDs by convolv- 
ing with the appropriate WISE filters, and we plot our 
synthetic photometry and redshift tracks in Figure El 

Colors combining the first three channels, 3.4, 4.6, 
and 12 /tm, are capable of selecting hyperluminous 
L m > 10 13 Lq ) gala x ies, particu l arly A GN/QSOs 
" 2011 IStern et all [20121: lYan et all 
'20121) found that the WISE color 



Eisenhardt et all 
20121 IStern et al 



cut [3.4] (ab) — [4.6] (ab) > 0.8 separates luminous AGN, 
which hav e been classi fi ed as such based on meeting the 
criteria in IStern et all (|2005l see Fig. EJ. As a com- 
plement to these techniques, we would like to separate 
luminous SF galaxies at high redshift. We find that only 
combining the first three WISE channels does not effec- 
tively separate our SF sources. The strongest separation 
of our tracks is produced by combining all four WISE 
channels as shown in Fig. [6l We overplot the AGN crite- 
rion of [3.4] (ab) — [4.6] (ab) > 0-8 as the grey dashed line. 
Our AGN redshift track confirms the IStern et al.l (|2012l ) 
criterion for AGN with z ~ 0.5 — 1.0 and z ~ 3. How- 
ever, none of our sources have a power-law slope steep 
enough in the appropriate wavelength range to meet this 
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criterion. 

Fig. [6] illustrates that the strongest separation of our 
sources and redshift tracks is for SF galaxies at z ~ 2 
due to prominent PAH features lying in the 22 /im trans- 
mission filter. The sensitivity dept hs at 12 and 22 
are ~ 1 and ~ 6mJy, respectively {Wright ct aL 2010). 
Given the sensitivity limits, at a redshift of ~ 2, a SF 
galaxy would need to be at least a ULIRG and proba- 
bly even a hyper-LIRG (Lir > 10 13 L Q ) to be detected 
by WISE. A benefit of this color selection technique is 
that few, if any, extragalactic sources with luminosities 
less than the ULIRG threshold should occupy the same 
region as the z ~ 2 ULIRGs. Therefore, our WISE color 
selection method is useful for selecting the brightest SF 
galaxies at z ~ 2. In the absence of WISE photome- 
try, Spitzer photometry can be used in similar combi- 
natio ns, with either 8 /im or 16 /im subs tituting for S12 
(e.g., Ilvison et al.|[200l IPope et al.|[200l . 

4.5. X-ray emission in AGN 

One of the motivating reasons for IR diagnostics of 
a galaxy's power source is to search for moderately to 
heavily obscured AGN missed by X-ray surveys. We use 
Chandra survey s of the GOODS fields (|Alexander et al.l 
I2003t ILuo et aLll2008D to determine if our diagnostics se- 
lect galaxies that are obscured and lack an X-ray detec- 
tion. For S250/S24 vs. Ss/S3„6, the X-ray detection frac- 
tion is 75% for AGN lying above the solid line and drops 
to 36% for AGN lying in the region dominated by SF 
galaxies. For Sioo/SW vs. Ss/S^.q, the X-ray detection 
fraction for AGN above the solid line is 79%, comparable 
to the detection fraction for the S250/S2A vs. Ss/S^.q di- 
agnostic. However, in the SF region, the X-ray detection 
fraction of AGN sources is 60%. Approximately 40% of 
the X-ray undetected AGN lie above the solid lines in 
both diagnostics, making our color selection criteria use- 
ful for selecting galaxies dominated by an AGN in the IR 
regardless of the level of obscuration. It is interesting to 
note that ~ 60% of our mid-IR dominated AGN would 
not be classified as such based on either X-ray data or 
the S250/S24 and Ss/S^.q colors. The AGN embedded 
in these galaxies are not having a strong enough effect 
on the far-IR and near-IR portions of the spectrum to 
distinguish them from SF galaxies on the basis of col- 
ors, although the AGN are having an effect on the mid- 
IR emission. This could be a product of viewing angle, 
or possibly even an evolutionary sequence where as the 
AGN grows more luminous, it will be reflected in its IR 
colors (see Paper I for a full discussion). 

5. APPLICATION OF NEW COLOR-COLOR DIAGNOSTIC 
TO ALL GOODS- HERSCHEL GALAXIES 

We now apply our new diagnostics defined above to 
broadly separate SF and AGN d ominated galaxies in the 
whole GOODS- Herschel survey (lElbaz et alJl20ll . We 
plot all galaxies in GOODS-N and GOODS-S detected 
in all four bands on our two new color-color plots in 
Figure [3 We use the color cuts derived in Sections 14.11 
and !4.2l (equations (3) and (6)) to separate sources dom- 
inated by AGN and SF activity in the mid-IR. Using the 
S250/S2A vs. Sg/S^.e plot (top panels), we have a total 
of 665 galaxies with detections in all four bandwidths in 
GOODS-N and GOODS-S, of which 58 (9%) lie in the 
AGN region. For Sioo/SW vs. Ss/S^.q (bottom panels), 



we have 988 sources with detections in all four badwidths, 
of which 94 (10%) lie in the AGN region. Both diagnos- 
tic plots lead to a similar fraction of GOODS- Herschel 
sources being AGN dominated in the mid-IR. 

It has been found that bright 24 /im flux density cor- 
relates with mid-IR AGN indicators at high redshift 
and in local ULIRG s (jDesai et a,l.ll2007t iDev et alJl25ol 
iDonlev et al.l 120081 respectively). Given this trend, we 
look to quantify the fraction of galaxies that are AGN 
dominated in the mid-IR as a function of flux density 
using our new color-color plots. In Fig. we plot all 
GOODS- Herschel galaxies with different colored sym- 
bols depending on their flux density (see legend for each 
panel). Our diagnostic was determined using 24/iJy 
bright sources, and applied to our own sample, the di- 
agnostic primarily recovered AGN with 6*24 > 300 /iJy. 
Furthermore, the contribution of the AGN to the mid- 
IR has been shown to increase with increasing 24 /im 
flux density, part icularly at higher redshift (z > 0.6 
iBrand et al.l 120061) . Therefore, it is perhaps not surpris- 
ing that only one fainter source (S24 < 200 /iJy) is found 
in the AGN region in Fig. 7. It is clear that a larger 
fraction of sources are in the AGN region of the color- 
color plot for the brightest 24 and 8/im flux densi- 
ties, whereas the 250 /im and 100 /im flux densities do 
not make as much of a difference to the AGN fraction. 
We quantify the fraction of sources in each flux bin that 
are classified as AGN and SF dominated based on these 
color-color plots in Tables [3] and HI Our results show that 
the presence of an AGN is only visible as an excess in the 
mid-IR, while the far-IR photometry alone is largely in- 
sensitive to the physics of the nuclear power source. We 
therefore conclude that the flux limit of a mid-IR survey 
has a big effect on the fraction of sources which have a 
significant AGN whereas far-IR/submm fluxes correlate 
less strongly with AGN fraction. 

Using the data in Table [31 we plot the percentage of 
galaxies which are classified as AGN and SF greater than 
a given 24 /im flux density in Figure [8] The fraction of 
AGN sources shows a steady increase with 24 /xm flux 
density. Above S24 ~ 750 /iJy, the majority of galaxies 
selected at these wavelengths will be AGN dominated. 
In Table [TJ we compare the intrinsic luminosities of our 
composite SEDs. Despite our AGN SEDs having an Lir 
below that of the z ~ 2 SF SED, the AGN subsamples 
have a higher median S24 flux density, particularly the 
featureless AGN. This effect is due to the differing ratios 
of far-IR to mid-IR luminosity. With the possible excep- 
tion of z ~ 2, where the brightest PAH complex falls in 
the 24 /im filter, at all redshifts AGN should comprise the 
majority of galaxies above S24 > 750 /im despite possibly 
being less intrinsically luminous in the IR than SF galax- 
ies (see Paper I for a direct comparison of our SEDs). 

6. CONCLUSIONS 

We have combined deep photometry from 3.6 — 500 /im 
with Spitzer mid-IR spectroscopy for 151 high redshift 
(z > 0.5) (U)LIRGs in order to explore the relative po- 
sitions of high redshift mid-IR AGN and SF galaxies in 
IR colorspace. We start by applying color-color diag- 
nostics based on IRAC colors alone (|Lacv et al.l 120041 : 
iStern et all 120051 ) to our high redshift sample and find 
that there is significant contamination of SF sources in 
the AGN regions; however, our mid-IR spectroscopic 
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Fig. 7. — Our new color diagnostic for separating AGN and SF galaxies presented in Figs. [4]and[5] (top and bottom panels, respectively) 
applied to our full sample of MIPS 24 fim selected galaxies in the GOODS fields which have corresponding Herschel photometry taken in 
the GOODS- Herschel survey. In the left panels, we color code galaxies according to their 250 /xm flux density (or 100 /im flux density); 
in the middle panels, according to 24 /im flux density; and in the right panels, according to 8 fim flux density. The region below the solid 
lines are the same as we defined in Figs. [4] and [5] a source lying above the line we label as AGN-dominated in the mid-IR (based on our 
extensive mid-IR spectroscopic sample). Our new diagnostics indicate that roughly ~8% of the GOODS- Herschel galaxies have an AGN 
dominating the mid-IR luminosity. Breaking things down as a function of flux density, we find that sources with the brightest 524 and Sg 
are predominantly AGN, whereas 5250 an d 5ioo do not preferentially select AGN at brighter flux densities (see Tables \S\ and [4j . 



TABLE 3 

Percentage of AGN and SF field galaxies according to S250, S24, and Ss strength. 





0350 






S24^im 










Range (mjy) 


% AGN (#) 


% SF (#) 


Range QuJy) 


% AGN (#) 


% SF (#) 


Range QuJy) 


% AGN (#) 


% SF (# ) 


< 15 


5 (25) 


95 (480) 


< 200 


1 (3) 


99 (352) 


< 40 


1 (6) 


99 (497) 


15-30 


12 (21) 


88 (150) 


200-600 


6 (19) 


94 (278) 


40-200 


10 (17) 


90 (151) 


> 30 


26 (12) 


74 (35) 


> 600 


51 (36) 


49 (35) 


> 200 


67 (35) 


33 (17) 



Note. — Data corresponds to the top panels of Fig. [7] which shows GOODS-N and GOODS-S field galaxies plotted on our color-color 
diagnostic defined in Fig. [4] AGN are those lying above the solid line indicated on the plot. In this table, we calculate the percentage 
(number) of AGN and SF sources in a given flux density range for S25O! S24, and Ss- 



sample confirms th e new IRAC color diagnostics from 
Don lev et all pOll . Adding 24 fim photometry does not 
effectively separate all of the AGN, but it does produce a 
separation of the SF galaxies according to redshift. The 



same separation of SF galaxies is produced by plotting 
our redshifted SED tracks in color-color graph using the 
WISE transmission filters. 

Our high redshift AGN and SF sources exhibit different 
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TABLE 4 

Percentage of AGN and SF field galaxies according to Sioo, S24, and Ss strength. 



Range (mjy) 


% AGN (#) 


% SF (#) 


Range (fjjy) 


% AGN (#) 


% SF (#) 


Range (^tjy) 


% AGN (#) 


% SF (# ) 


< 5 


5 (40) 


95 (779) 


< 200 


3 (22) 


97 (613) 


< 40 


3 (21) 


97 (741) 


5-15 


15 (29) 


85 (165) 


200-600 


9 (35) 


91 (337) 


40-200 


13 (33) 


87 (222) 


> 15 


36 (25) 


64 (44) 


> 600 


49 (37) 


51 (38) 


> 200 


62 (40) 


38 (25) 



Note. — Data corresponds to the bottom panels of Fig. [7] which shows GOODS-N and GOODS-S field galaxies plotted on our color 
diagnostic presented in Fig. [5] AGN sources (94 total) are those lying above the solid line indicated on the plot. In this table, we calculate 
the percentage (number) of AGN and SF sources in a given flux density range for Sioth S24, and Ss. 



spectral features and SED shape in the near-, mid-, and 
far-IR/submm wavelengths. We explore combining pho- 
tometry from all three IR ranges for an optimal selection 
technique. In addition to photometry, we explore where 
AGN or SF sources should lie on average in different color 
combinations by using redshift tracks of our composite 
SEDs. We present two new color-color diagnostics com- 
bining Spitzer mid-IR and Herschel far-IR photometry 
that can be used to estimate whether a galaxy harbors 
an AGN in the absence of IRS spectroscopy. The opti- 
mal color diagnostics for bright S24 sources spanning the 
full IR spectrum are S250/S24 vs. Ss/S^.e and Sioo/SW 
vs. Ss/Ss.q. These diagnostics have a low contamina- 
tion rate of ^10% in each of the SF and AGN occu- 
pied regions, and the contaminating AGN chiefly have 
S24 < 300 /iJy. Moreover, the S250 /S24 color estimates 
if an AGN is significantly contributing to the bolometric 
output of a galaxy. The mid-IR photometry from WISE 
can also be used in conjunction with Herschel photom- 
etry to separate out AGN. When only limited data is 
available, either of the colors S250/S24 or Ss/S^.e can be 
used alone to select AGN. 

We apply our new color-color diagnostics to the entire 
GOODS- Herschel survey to determine the fraction of 
sources which are dominated by AGN and SF activity. 
We find roughly 10% of GOODS- Herschel galaxies have 
a significant AGN. We also confirm that a higher frac- 
tion (> 50%) of the brightest sources, Ss > 200 [iJy or 
S*24 > 600 jitJy, have colors indicative of a bolometrically 
significant AGN. The 100 /im and 250 /im fluxes show a 
much weaker dependence on AGN fraction, though when 
combined with S24, 5*250 effectively separates out the 
AGN, indicating that the amount of cold dust emission 
relative to mid-IR emission is a good indicator of the IR 
power source. We conclude then that far-IR fluxes or 
colors alone cannot be used to determine the nature of 
the IR power source, and the most effective method is 
to combine mid-IR and far-IR data. 



100 - 




S 24 (nJy) 

Fig. 8. — We calculate the percentage of AGN and SF sources 
(with Poisson error bars) greater than a given 24 fim flux density. 
Sources are taken from the GOODS fields and defined as an SF or 
AGN galaxy based on their position in S250/S24 vs. Sg/Ss.e color 
space (see top panels of Fig. 0. The brightest flux densities are 
dominated by AGN, illustrating that applying a flux density cut 
in S24 is useful for selecting a population with a large percentage 
of AGN. 
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